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To study composite fermions around an even denominator fraction, we adapt the phase space 
£N| , single-particle path integral technique for interacting electrons in zero magnetic field developed 

recently by D.S. Golubev and A.D. Zaikin, Phys. Rev. B 59, 9195 (1999). This path integral 
description gives an intuitive picture of composite fermion propagation very similar to the Caldeira- 
Leggett treatment of a particle interacting with an external environment. We use the new description 
' to explain the origin of the famous cancellation between the self-energy and the vertex corrections 

in semi-classical transport measurements. The effective range of the cancellation is given by the 
size of the propagating particle, which for the Coulomb interaction scales with the temperature T 
as T _1 / 4 | logT| _1 in the semi-classical limit. Using this scheme we find that the effective mass in 
the semi-classical limit for composite fermions in GaAs is approximately 6 times the bare mass. 
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I. INTRODUCTION 



In this aeticle we address the question of semi-classical propagation of an interacting system in a real space 
formulationEI. In particular we have in mind the composite fermion (CF) system which is seen in the fractional 
quantum Hall regime around filling factor |. The experimental efforts so far indicates that this system is to a high 
degree a classical system in-.the sense that it is very hard to observe any kind of interference effects. 

In elementary textbookso the notion of semi-classical transport is sometimes introduced through localized wave 
packets, and these appear to be essential building blocks. Historically this has really not always been the case. For 
instance the derivation of the Boltzmann equation in the case ofrimpurity scattering typically only uses individual 
quasi-particles with a given momentum i.e. delocalized plane wavesu. This is also the case for other weakly interacting 
t-H , systems like electron scattering on phonons at low temperatures. However, a Boltzmann equation for electrons 
■ scattering on hot phonons could not be derived using plane waves. Prange and Kadanoff solved that problem by 
shifting focus from individual quasi-particles to deformations of the Fermi surfacecl. In real space this step corresponds 
to forming an electron wave packet that is localized in the direction of propagation. This localization in the direction 
of propagation causes the electron-phonon scattering to be very brief in time due to the huge difference between 
the Fermi velocity and the sound velocity. Recently Kim, Lee, and Wen showed that the Boltzmann equation for 
' composite fermions only makes sense for smooth Fermi surface deformationsQ. In real space the wave packet now 
f^i \ also needs to be localized in the direction transverse to the direction of propagation, so the composite fermions is the 
"j""^ ■ first system where the textbook introduction actually is a proper description. For composite fermions the necessity 
of transverse confinement of the wave packet can be understood intuitively. The important interaction is mediated 
through a gauge field and the coupling is thus through the flux enclosed by an area. For a localized wave packet with 
a finite width propagating along a classical path the relevant area is the area swept by the wave packet. The width 
of the packet is thus a natural cut-off on the transverse gauge fluctuations. In this article we pursue some of the 
consequences of this idea. 



, We will study the real space propagation with the help of phase space path integrals. Originally, the Fevjiman path 
L* ■ integral technique was developed for non-interacting systemsa. Feynman-VernonFj and Caldeira-LeggettQ expanded 
the idea to systems interacting with an environment. Recently Golubev and Zaikin extended the technique to cover 
the linear response regime for the reduced one-particle density matrix in a Coulomb interacting gas of electrons!!. 
This last step is not yet widely recognized. Nonetheless, in Sec. [jl] we ad apt i t to the case of composite fermions. We 



summarize the main points of the path integral technique in subsection II G. In Sec. [II the path integral technique 
is on a formal level compared with the Kubo formula description of transport and, more briefly, with the Landauer- 
Buttiker scattering approach to transport. In Sec. |y| we discuss the typical semi-classical approximation scheme of 
saying that the two paths are almost identical, and the typical width between the two paths is calculated. In Sec. |v| 
the calculation of the effective mass of composite fermions is carried through in the semi-classical limit. 



*Present address: Nordita, Blegdamsvej 17, DK-2100 Copenhagen, Denmark 
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We have found it necessary to introduce quite a few averages and short hand notations. As a readers guide we list 
them here. (X) is the full quantum and thermal average. A subscript as in (X) means that the average is performed 
for the non-interacting electron gas. (X) a t is the average over different Chern-Simon gauge field configurations, the 
precise definition is Eq. (|l8|). \X) = \X)\ a in Eq. ( |78| ) is the average in the time interval [0,t] over different 
single-particle paths, where each path is given a weight that reflects its importance in a given experiment. It thus 
depends on what experiment is considered. (X) = {X)ip,t] is the shorthand notation for the sum over single-particle 
paths including the weight due to the interactions. It is defined in Eq. j7S|). ((AT) defined in Eq. (A4) is a short hand 
notation for the Grassmann integral. It is only used in Appendix |A|. Quite generally the superscript 'eq' refers to the 



average being calculated in equilibrium. For instance (A)° q t in Eq. (pll) is the equilibrium version of (X) a t . 



II. THE DENSITY MATRIX 



We consider a two dimensional electron gas in a static and homogeneous external magnetic field B cxt perpendicular 
to the plane. We write B cxt as a scalar since its direction is fixed, and likewise we consider the rotation of any 
vector potential in the two dimensional plane as a scalar. Furthermore, we also apply an external electrical potential 
4>(r,t). The external potentials are written in four-vector form as A — ((f), A). In this section the goal is to derive a 
path-integral description of the linear response of the reduced one-particle density matrix p to the potential 4>(r,t). 
For simplicity we will in this section neglect the influence of impurities, but it is straight forward to include them in 
the formalism as an extra external potential. 



A. The formulation of the problem using electrons as the fundamental particles 

We will assume the system is described by the Hamiltonian 

H cl = H a [A] + H int . (1) 

Here 

- e4>(t,t) - h^(t) (2) 

is the Hamiltonian for non-interacting particles in the 2D plane represented by the second quantized electron field 
operators ip(r) and V^(r), while 

#int = y fdrf dv'^(v)^(v')v(v - r')$(r')i>(r) + H ch (3) 

is the part due to electron-electron interaction. 7J e b represents the interaction between the electrons and a fixed 
neutralizing positively charged background jellium. m is the mass of the electron, /i is the chemical potential, — e 
is the electron charge, and e 2 v(v) is the Coulomb pair interaction between electrons. If the screening effect of the 
back gate is neglected v(r) = l/(47re|r|), where e is the permittivity of the medium surrounding the two dimensional 
electron gas. 

All the physical one-particle quantities of the system, like the current and the density, can be expressed in terms of 
the reduced one particle density matrix 

p c (r,r';t) = (ft(r',t)i>(r,t)) (4) 

Here (X) denotes the full quantum mechanical and thermal average of X. The subscript c stands for canonical, 
because upon a Fourier transformation the annihilation and creation operators correspond to states labeled by the 
canonical momentum. This labeling is not gauge invariant so neither is p c (r, r'), but it has some nice formal properties 
we will make use of in the following. Final physical expectation values we prefer to phrase in terms of the gauge 
invariant density matrix 

p(r,r')=Pc(r,r')e i *X' dr "- A(r " ) ! (5) 

where the path of integration is chosen to be along the straight line from r' to r. In general different choices of the 
integration path lead to different values of p, but we are mainly interested in local probes like the density and the 
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FIG. 1. The Grassmann fields i/S e i, ip e i, ■ip, and tp and the Chern-Simon gauge field a = (ao,a) are all live on the contour Ct- 
As can be seen from the drawing the contour Ct runs along the physical time axis from — oo up to time t and then back again 
to — oo. The contour time a is on the forward and backward part of the contour denoted a = (s, +) and a = (s, bi) respectively, 
where s is the physical time. To each physical time s < t there is thus as an example two independent Chern-Simon gauge 
fields a(s,fin) and a(s,bi). If in particular t = oo the usual Keldysh contour Coo is recovered. 



current and then any path that locally approaches a straight line yields the same result. In such local probes it is 
convenient to introduce variables for the average R = |(r + r'), and for the difference Ar = r — r'. We use a new 
symbol for the density matrix expressed in these variables 

/ Ar Ar\ 
u ,(R,Ar) = p(R+— ,R-— J- (6) 

The Fourier transform with respect to the difference variable is the corresponding gauge invariant Wigner distribution 
function 

w(R,p) = J dAr e~^ p ' Ar w(R,Ar). (7) 
In terms of the gauge invariant Wigner distribution function the current density at R is given by 

Notice the momentum in the Wigner distribution p in the classical local limit is the kinetic momentum. 

The thermal and quantum average of p c and p can be calculated with the help of a path integral over Grassmann 
fields and ip \. -0 c i corresponds to the operator ^ and they should not be confused. Likewise, ip e i corresponds to 
if). Using this representation we obtain 

( '-t)= 'f f, rfr"-A(r") / P^cl / X>^cl^el(r / , t)lp cl (r, t) exp (j-S[{p ch V>ei; tj) 
Mr ' r ' j_e /X^el/^eiexp^^Vei;*]) ' 

where the action S is given by 

S[ipei, ip e i;t] = J da J dr ip c i(r, a)ih-^-ip c \{v , a) - J daH c \[^ ch ip c \\. (10) 

The contour Ct runs from minus infinity up to the observation time t and back again to minus infinity. The contour 
time a = (s, /x), where s is a real time and /i = ± is a contour index that takes the value + on the forward part of the 
contour and — on the backward in time part of the contour, see Fig. [|. In the prefactor ■iA e i(r / , t)-0 c i(r, t) the creation 
Grassmann field ip e \(r',t) has to come later than the annihilation Grassmann field ip c \{r,t) on the contour Ct, but 
whether both fields are on the forward or backward parts of the contour Ct or ip e \(r,t) is on the forward part, while 
tp e i(r',t) is on the backward in time part of Ct is immaterial. 



B. The Chern-Simon field dependent density matrix 



The first step towards obtaining an effective single-particle path integral description is to introduce fluctuating fields 
as extra integration variables such that the action becomes quadratic in the Grassmann fields. The electron physics is 
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then reduced to that of non-interacting electrons in a fluctuating field background. In the case of zero magnetic field 
the relevant transformation is a Hubbard-Stratonovich transformation that introduces an effective electric potential at 
each contour pointtl. In the case of composite fermions the relevant transformation is a Chern-Simons transformation 
that introduces a Chern-Simon gauge field a — (ao> a ) at each contour point, and that also performs a singular gauge 
transformation on the Grassmann fields attaching two magnetic flux quanta h/e to each electron. The transformed 
Grassmann fields we denote by %p and ip, not to be confused with yji, ib e \ and tp e \. 
The action transforms into the standard action for composite fermionsu 

Scf [i>i V*) a > Arff — o,; t] = Scs [o; t] + e f da f drao(r,a)n + j da j drip(r, a)ih—?p(r, a) — f daHa [A c g — a] . 

JC't J Jc t J "° Jc t 

(11) 

Here jSosfflji] is the Chern-Simon action 

f f ( e 2 e 2 d&\ 

S CS [a;t] = J daj dr [^a (a)V x a(a) + — a(<r) x - j 
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J da J dr J dr' [V x a(r, a)} v(r - r') [V x a(r', a) 



(12) 



The average mean- field magnetic field 2nh/e has been subtracted from both the externally magnetic field V x A 
and the Chern-Simons gauge field. The remaining fields are denoted by A e g — (0, A e ff) = (0, A — (Vx) 2— n) and 
a = (ao , a) respectively. For notational consistency with the subsection on linear response we have chosen to have 
a vanishing external potential, = 0. If non-zero it would be added to A e g. The last term in Scs originates from 
the interaction term H mt , Eq. (||), after applying the identity V x a(r, a) = 2^ ( , tp e i(r,a)ip e i(r,a) — n). For later use 
we notice that in the terms that couple the Chern-Simons gauge field a to the Chern-Simons fields tf> and -0, a only 
enters in the combination A Q g — a. 

In the first line of Eq. ( |i~l|) the temporal integration contour can be extended to the Keldysh contour Coo, since 
without a Hamiltonian there is no real time development and consequently the integral over the a-field on the Keldysh 
contour from the observation time t to infinity is one. is used below for technical reasons whenever the evolution 
of various quantities like effective actions and density matrices in a fixed background gauge field history a is discussed. 
The contour Ct is used when the connection to a physical density matrix as in Eq. ( |l5|) is wanted. 

The prefactor in Eq. (|^) transforms as 

M* f ,t)M*,t)e i *£* r ''' M ''' ) ^ ^(r , J *)^(r,t)e i *^*"-^ (, "' t) - a(r "'* ) ] J (13) 

where 

a(r",t) = i[a(r" ) t,+)+a(r",t,-)]. (14) 

is the average of the Chern-Simons gauge field on the two parts of Ct- This average value is chosen for notational 
consistency. It is allowed since the a-field is continuous on the contour Ct, so it does not matter whether the Chern- 
Simons field at the turning point t is evaluated on the forward or on the backward part of the contour. The total 
expression for the density matrix in terms of the transformed fields is 



p(r,r';t) 



i-ScF 



Va 



V{p f PV^(r',i)^(r,t)e^ ScF+lf /> r "-[ Arff(r " Kii(r "'* ) L (15) 



The exponent is quadratic in the Grassmann CF variables ip and ip and thus describes non-interacting particles in 
a gauge field history A e g — a. A typical gauge field history is not physically realizable since the gauge field a at a 
given time s is different on the forward and backwards parts of the contour Ct, a(s,+) =/= ct(s, — ). Despite this we 
will in analogy with the usual definition of a density matrix Eq. (Q) and following Golubev and ZaikinH define the 
Chern-Simon field dependent density matrix 



Pac{A c ft- a;r,r',t) 



Pr I Pi'vxp ( ^Scf 



ZV / Vi; i/>(r>,t)il>(T,t)exp [ ^S CF ) . ! U>) 



Notice that the right-hand side only depends on the gauge fields in the combination A c g — a, since the contribution 
from 5cs[a;^] to the action Scf cancels out. The subscript a on p ac is to remind us that p ac depends on the field 
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history, while the subscript c stands for canonical in analogy with p c , see Eq. (J4J) . The physical density matrix is 
recovered from p a c(A e s — a; r, r', t) by averaging over the different a-field histories 

p(r,r',t) = ^p ac (A eS -a;r,r',t)ex.p^£ i dr" ■ [A eff (r") - a(r", t)]j^ . (17) 

In the average 

_ JVaXe X p(±S oB [a;t]) 



the weight 



exp ( -S e s[a;t) j = T>if> / Vipexp ( -Sc F [^,<0,a, Aiff - a;t] ) . (!.!)) 



is the denominator in Eq. (|i~6|). The effective action S^fa;^ for the a field also depends on the external gauge field 
A e s, but this dependence has been suppressed in the notation, since it is not needed in the following. Eq. ( |17| ) is an 
important formula, since it establishes the connection between the density matrix p ac depending on a given history 
of the gauge field and the physically accessible density matrix p. 



C. The effective action for the a-field 



Introducing the a-dependent current j a (^4 c ff — r, s) at time s 
j a (A, ff -a;r, s 



—e 
2m 



-^- + eA off (r) - ea(r) - + eA cff (r') - ea(r') 

i Or i or' 



p ac (r,r';s) 



it is shown in appendix ^ that the effective single-particle action S c g [a; t] is given by 

S e tt[a;t] = Scs[a;t] - / ds dr [eAa (r, s)(p ac (r, r; s) - n) + Aa(r, s) ■ j (r, s)} 

J — oo J 

where 

Aa Q (r, s) = a Q (r, s, +) - a a (r, s, -) 



(20) 



(21) 



(22) 



is the difference of the Chcrn-Simons gauge field on the two parts of Cj. 5cs[a>*] is the usual Chern-Simons action, 
see Eq. (|l2|). In the notation the explicit dependence of the density matrix p ac and the current j a on the gauge field 
A e f{ — a has been suppressed. In equilibrium, which is denoted by the superscript " eq" , we have to second order in 
the a-field 

S^[a;t] « S cs [a;t] - £ /* ds' f ds" f dr' f dr"Aa a (r', s')R a0 (r', s'; r" , s")a (r", s") 

a 3 J —oo J J 

ds ' f ds" J dr' J dr'Aa a (r', S ')F a0 (r' lS ';r",s")Aap(r", S "). (23) 



a , ' 



The summation variables a, (3 £ {0, x, y} run over the single time-like and the two spatial indices e.g. a a = (ao, a x , a y ). 
The terms linear in a vanishes, since p ac (A cB ,r,r,s) = n and j a (A c ff; r, s) = 0. The retarded 3x3 response matrix 

(^Rj is given by the usual Kubo formula apart from an extra minus in the part that mixes the density with the 



currents. 



SfV «'V ,"1 = _1 ( {[~ep ac {r' , s'),-ep ac (r" : s")}) - {[-ep ac (r', s'), j a (r", s")]) Q \ , ,_ 
1 ' ' ' n\ -([ja(r', S '),-ePac(r", S ")]) ([ja(r', s% j a (r", s")}) ) U[S S > 



+ 











<L- Pac (r',s')S(r'-r")6(s'-s") 



(24) 
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and the so-called fluctuation kernel 



(f) by 

1 X ' 11 2 ' 2j ~ 27i ^ - ({ja(r', *'), -ep ac (r", s")}) <{j a (r', S '), j a (r", S ")»o J ' W 

The dependence of the density matrix p ac and the current j a on the gauge field A c g — a has been suppressed in the 
notation. (X) Q means the average over the non interacting electron gas without the Chern-Simons gauge field a, and 

is the step function. In equilibrium the retarded response matrix R is related to the fluctuation matrix F through 
the fluctuation dissipation theorem 

F al3 (r i , r 2 ; w) = l - coth ( — — ) [R a/3 (n , r 2 ; u) - Rp a (r 2 , ri ; - ui)} , (26) 



2 \2k B T 

where u is the Fourier frequency corresponding to the time difference s' — s" . It is used that in equilibrium R and F 
only depends on s' — s". In terms of the 2x2 conductivity tensor <r(q, uj) the response matrix R is written as 

% )W ) = if^ q :f( q 't ;) ' q - q :f (q \ w) V (27) 

V ~cr(q,w)-q w<r(q,w) / 
Below we are mainly going to need the retarded gauge field propagator D a p defined by 



D ap (r', a'; r" , s") = --e 2 (a Q (r', s') A a/3 (r", S "))° Q q t = --6 (s' - s") e 2 <[a Q (r', S ') , a/3 (r", s ")]) (28) 



and the correlation function C a p for fluctuations of the field 

C a0 (r', s'; r", s") ee ^ (a a (r', S ') ^ (r", S ")) C a q t = ^~ ({a a (r', s') , a (r", «")}>„ , (29) 
while the auto-correlation function for Aa a vanishes: 

(Aa a (r',s')Aap(r" lS ")) c a l = Q. (30) 

The average 

_ Jl? a Xexp(^ c >;t]) 



is the average (X) t of Eq. (18), calculated in equilibrium. 

In the Coulomb gauge the most troublesome cpnttjibutions arises from the low frequency limit of the so-called 
transverse part of the gauge field propagator Z)^ Ji3'E!J, which in Fourier space is given by 

£ ±± (q^) = ^V 7 *§^(q, W ). (32) 
ijkl V 

Here e xx = e yv = 0, e xy = —e yx = 1, and the indices i, j, k and I run-over 'x' and 'y'. In equilibrium in the limit 
hw <C qvF the propagator D±± is given by the standard approximation!!!! 



D L± ( q. ) 2 :rh±- 9 ! a „,„, n 3 = ^TTh j- j ] . (33) 



PF iuJ _ 2 */v( q )q* p F iuJ _ g q r, 

For the Coulomb interaction v (q) = ^ so r\ = 1 and gi = ^"t^V^c m terms of the average Coulomb energ; 



?y 



ec = j^V™- n is the density of the electron gas. Other values of 1 < rj < 2 corresponds to a power law decay r n . 
The fluctuation correlation function C is given by 



C(l,2) = j d3 J d4 5(1,3) F (3,4) D (2,4) , (34) 



G 



where 1,2,3 and 4 each are a short hand notation for a full set of variables a,r,s. In equilibrium C is determined from 
the fluctuation dissipation theorem ( |26| ) 

C^(r',r";o;) = |coth(^) [£> Q/3 (r',r'» - D^ a (r",r'; - w )] (35) 

In particular in the Coulomb gauge is the most singular contribution due to the transverse part of the propagator 



f „ , ( q. ^ i % | 6 a0 - ^f-J coth [^-^ ) rm/> i i q. - i I 30 i 



D. The equation of motion for the gauge field dependent density matrix 

In appendix |A| we have by brute force differentiation showed that p ac fulfills the all important nonlinear equation 
of motion 

iftdtPac = (1 ~ Pac)H(t, +)p ac - p ac H(t, -)(1 - p ac ). (37) 

This equation should be read as an operator identity, i.e. we have suppressed all arguments and the spatial integrals 
over the internal variables, 1 is a shorthand notation for the identity operator S(r' — r'"), H(+) is the kernel in the 
quadratic action for the Grassmann CF variables on the forward part of the contour, i.e. 



H(v', r"; t, +) = 8(r' - r") J -e<j>(r",t) + ea (r",t, +) + ±- 



-V" + eA cff (r")-ea(r",i,+) 

i 



(38) 



a(+) = [ao(+),a(+)] is the gauge field on the forward contour. Likewise is H{—) the kernel on the backward in 
time part of the contour. If the gauge field is identical on the forward and backward sections of the contour Eq. 
( |37| ) specializes to the usual equation of motion for the density matrix iTidtp = [H, p}. The extra complication in the 
general case is the term quadratic in the density matrix: —p ac [H(+) — H(—)\p ac . However we restrict ourselves to 
the linear response regime. 

E. Formulation of linear response 

We are interested in the response of the system to an external driving field </> and in particular in the change in 
the density matrix: <5p(r, r';t) = p(r,r';t) — p cq (r, r';t). The superscript "eq" still denotes equilibrium where </> = 0. 
In general both the gauge field dependent density matrix p ac {A e ft — a; r, r', t) and the effective action for the Chern- 
Simons field S e f{[a', t] depends on 0, giving rise to two different contributions to Sp(r, r', t) — compare with Eqs. ( [T^ ) 
and (|l^). The contribution coming from the variation of S c g[a; t] can be accounted for by a change of variable from 
the externally applied field 4> to the total gauge field felt be the composite fermions 

AtotO; t) = {4>{t),0) + A ind {a; t). (39) 

Here Aj n( j is the self-consistently induced gauge field due to the change in the current and the density when 4> is 
applied. Technically it is defined by the demand that 

S eS [a- A ind (a);i\ = S°%[a;i\. (40) 

From this definition it is not immediately clear that the induced field is the same on the forward and backward part 
of the contour and thus can be interpreted as a physical field. However a careful inspection of Eqs. ([[2]) and (|2l| ) 
shows this is the case. Furthermore it is seen that the induced electromagnetic fields corresponding to Aj n( j are 

47T?l 

Bhid{A cS - a; r, t) = Sp ac (A cS - a; r, r; t) (41) 

e 

Eind^cff - a;r,t) = —^z x 5j a (A cS - a;r,t) + eV / dr' v(r - r')Sp ac (A cS - a;r,r;t), (42) 
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where 



Sp ac (A cS - a; r, r', t) = p ac (A eS + A tot - a; r, r', t)e^ £' dr " Atot - p c £(A cS - a; r, r', i) (43) 
Sj a {A cS - a; r, t) = j a (A off + A tot - a; r, t) - f a q (A cS - a; r, t) (44) 

are the changes in the density and the current respectively, induced by the total field A tot . The phase factor 
exp [if J , dr" ■ A to t] has been introduced in order to let Sp ac refer to the equilibrium state. For instance S_] a re- 
sults, when 5p ac is substituted for in the expression for the equilibrium current Eq. (|2p|). Eqs. ( [4l|) and (42) are 



what is expected on physical grounds. There is a magnetic field associated with a change in the density of composite 
fermions and thus the amount of attached flux. There is also a Faraday electric field associated with a current of 
attached flux. The last term in Eq. fl42[ ) is due to the electron-electron interaction. 

The Jacobian associated with the variable change from a to a — A- ln ^(a) in Eqs. jl7| ) and ( |l8| ) is the identity 1, since 
both the density and the current, and thus A- m d(a), are retarded functions of the gauge field A c g — a. Consequently 
the change in the physical density matrix 

5p(r,T';t) = ^ ac (A e ff-a,r,r';i)exp|i|^ dr" ■ [A eff (r") - a(r", t)] j ^> (45) 

is the average of the change in the a-dependent density matrix over the different a-field histories with the weight 
exp [fS'offH] pertaining to equilibrium. 

Golubev and Zaikinlil did not perform this change of variables, since the difference between the external and the 
total field is not of great bearing in their application of the formalism to weak localization. The second term in their 
formula (46) which they neglect is exactly what is covered by this change of variable. For composite fermions such a 
neglect is not permissible, since the extra term originating from the variation of S^s [a; t] gives rise to for instance the 
Hall conductance. In our formulation the Hall effect arises from the electric field induced by the current — the fipt 
term in Eq. (^) . The mechanism is completely analogous to the treatment in the case of the Boltzmann equatior£3. 

Notice that so far no linearization in the applied field has been performed and everything is exact. 

The linearized equation of motion for Sp ac (A c ff — a; r, r', t) is obtained by linearizing Eq. ( |37j ) in the total external 
field 

ihd t 8 Pac (v, r'; t) ~ J dr"H+{r, r", t)6p ac (v", r'; t) - J dr"6p ac (r, r"; t)H-(r", r', t) + ihD(r, r', t), (46) 

where the explicit dependence on the gauge fields A c g — a and A tot is suppressed in the notation. New one-particle 
Hamiltonians 



H+= H(+) + H(-) 



4=0 



\-Pac)[IH )-././M 



(47) 

6=0 



H~ = 



H(+) + H(-) 



[H(+)-H(-)} 



= 



Pat 



(48) 



0=0 



have been introduced in analogy with the equation of motion for the physical one-particle density matrix ifid t p = [H, p] . 
H + =/= H~ because the Chern-Simons field is different on the forward and backward contours. The inhomogeneous 
stimulus D = Di + D 2 is split into two parts, each of which are gauge invariant in the total applied field A tot . D 2 
contains all the dependence upon the difference gauge field Aa, while Di is independent of Aa. In Appendix |b] the 
explicit expressions for Di and D 2 are listed. In linear response Di and D 2 contribute additively to the physical 
density matrix. The response due to D 2 we will neglect, since it vanishes in the approximation where D is replaced 
by its average value (D)° q t . Below we do not go beyond this approximation. 

In the real space representation used above D is not gauge invariant in A c g(r) — a, so we find it convenient to change 
representation to the gauge invariant Wigner function set of variables. In analogy with Eqs. (||), (||) and ([?]) 



D(r,r',i) = 



dp 



■ J rl dr"-[p-eA off (r")+ea(r",t)]-p 



2^ + r '):P 



(49) 



When the generalized density matrix p ac (r,r') is peaked around r = r' on a length scale which is much smaller than 
the typical scale of variation of the total driving field, the gauge invariant Di and thus D is approximately equal to 
the standard Boltzmann equation driving term 
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D(R,p) 



1 \ d 

> [ — p X zBtot + Etot ) ' — W a (A eS - d] R, p). 



dp 



(50) 



In the presence of the gauge field A g — a the Wigner distribution function w a (A e s — a; R, p, t) is defined in analogy 
with Eqs. (|) and @) 



w a (A cS - a; R, p,t) = / dAre 



i In- A *"-[p-eA (>tf (r")+ea(r",i)] 



Pac(A cff - a;R+ -Ar,R- -Ar;t). (51) 



In the semi-classical limit it corresponds to a field dependent phase space distribution function. We will not consider 
effects beyond the approximation where it is the Fermi distribution function. 

F. The single-particle path integral in linear response 

In this subsection we obtain in terms of a single-particle path integral the formal solution of Sp in the linear response 
i.e. Eq. The first step is to notice that the solution to Eq. ( |4^ ) can be written in terms of evolution operators 



5 Pac (T,r';t) = ~j^J dt* J dr x J dr 2 [7+(r,t,n,t i )D(ri,r 2 ,t i )C7Z(r 2 ,t r',t). 



(52) 



The evolution operators 



U»(t,t') = f^exp 



h.L, 



dt"H' 1 (t"" 



fj, = ±, 



(53) 



where T + and T are the time and anti-time ordering operators respectively. 
The corresponding change in the physical Wigner function is 



£u;(R, p, t) 



dt; / g?R, 



dpi 



{2irhf 



(J(R, p, t; Ri, p^ fj)D(Rj, pi, i;)) ' 



eq 



(54) 



where 



-f r R / dr"-fp-aAeff(r")+ea(r",t)l I f * Ar- <ir"-[pi-eA cff (r")+ea(r",t i )l 



Ar Ar Ar 4 Ar» , 

xJ(R+— ,R-— ,i;R+^-,R i -^ ) t i ) 

is the propagator of the stimulus D in the Wigner representation. In the real space representation 



(55) 



1 7(R+-^,R-^,t;R+^p,R i -^£*,t i ) = cr+(R + ^ >t; R i + ^ )ti )[r(R i -^,t i ;R-^pf). (56) 



The evolution operators £/ + and U can be represented as usual single-particle phase space path integrals 



^(ri,ti;r 2 ,ta) 



r(ti)=n 



r(t 2 )=r 2 



Vr Vpexpi- / dt [(p - eAsff + ea) • r - F"(r(fl),p(«),s)] 



(57) 



where fi = ±. The Hamiltonian kernels are given by 



p 2 e 



tf+(r,p, s ) = ^ 
2m 



2 Aa 2 (r 



8m 



- - w a (r,p(s+),s) 



e0(r) + ea (r, s) 

f Aa - — • Aa(r, sf] + ih^^ ■ Vw Q (a; r, p(s+), 
V to / 2to 



H (r, p, s) = — 



p 2 e 2 Aa 2 (r, s) 
2m + 



8to 



Aao • Aa(r, s) 

m 



e<f>{r) + ea (r, s 
1 



w a (r,p(s+),.s) 



ift^^ ■ Vm a (r, p(s+), s). 
2 m 



(58) 



(59) 
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(Ri + ^Sp»,ii) 



a Q (r,P,A« = +) 



o(Rj,Pi,ii) (R/.p/.t/)!! 



(Rj - ^.Pjjij) a a (r,p, n = -) (R/ f^j P/j i/) 

FIG. 2. The closed contour path (r(cr),p(a)) starts out at a = (ti,+) at (Ri + Ap-,Pi) and runs along (r(s, +), p(s, +)) 
up to (R + ^,p) at contour time a — (t, +)• Here it runs along the straight line path from R + y to R - ^- with the 
constant momentum p and at constant a. From here it then runs back along [r(s, — ), p(s, — )] to (Ri — ^p-,pi) at contour 
time a = (ti, —). Finally it runs along a straight line path back to R; + with the constant momentum Pi and at a constant 
time U 



Here s + on both the forward and backward in time contour is a bit closer to the observation time t than s. The only 
place where this small time shift is important is in the Wigner distribution function w a (r, p(s + ), s), which is likely to 
be a rapidly varying function of the momentum p at the Fermi surface. Note also that p in the semi-classical limit is 
the kinetic momentum in the gauge field A e f{ — a. 

Using the path integral representation Eq. (^) for U + and U~ we finally arrive at the following contour path 
integral representation for the propagator J 

J(R,p,f,Ri,pi,ti) = J dArJ dAn J dr[ J Vr(+) J Vp(+) J Vr(-) J 2?p(-)exp (^S + ^AS) , (60) 



where 



So = <j> dar ■ {p(a) - eA cff [r(er), s}} - dsH Q [r(s, +), p(s, +), s] + jT dsH [r(s, -), p(s, -), s] (61) 

is the action for non-interacting particles in an external gauge field A e g. The closed contour is depicted in Fig. ^. All 
the contributions related to the Chern-Simon gauge field are gathered in the second term in the exponent of Eq. ( |60| ) 

AS = +e <j) {dr ■ a[r(a), s] — d<Tao{r(a), s]} 

+e [. ds {l~ Wa[r{s ' +) ' p(s+ ' +) ' s] } {^T^ ■ Aa[r(s ' +) ' s] ~ Aa °[ r ( s ' +)' s l 

+6 /. dS {l - w <>[r(s,-),p(s + ,-),s]J | P( '^" ) ■ Aa[r(s,-),s] - Aa [r(s, -), s] 
£ ds { Aa[r(,s, +), s} 2 - Aa[r(,s, -), s} 2 } 
Below we only consider the approximation obtained by factorizing Eq. (|54| ) 

MR,P,t) = J dti J dR, J ^ I (J(R,p,t;R,,p l ,^)) C a q t (D(R J ,p,,i l )>°a' 



(62) 



eq 
t • 



(63) 



Above we argued that the stimulus (D)° q in the semi-classical limit is the standard Boltzmann driving term, Eq. (|5C 
We will calculate the average of the propagator J perturbatively to second order in the Chern-Simons gauge field a. 
This procedure is gauge invariant. First notice that (AS)° q t — 0, since (a = 0)° q = and all averages with a Aa at 
the latest time vanishes. So, to second order in the Chern-Simons field fluctuations is 

<J(R,p,t;R J , Pl ,t t )) G a q t - J dAv J dAv, J dv[ J Vv(+) J Vp(+) J Pr(-) J PpRexp Us + ~S lnt ) , (64) 



where 



s . mt =^(Asr: t . (65) 
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To the same order it is permissible in the expression for AS to use the approximation w a (a; R, p, t) = 
w a {a; R, p, f)| a= n = w (p) = /( e p)i where /(e) is the Fermi distribution function. The contribution from the gauge 
field to the effective action is then 



27? 

i 
2 



fi'n" 



ds' I ds" /"V ( -1 , * (*') M') ) C (r 0', m') , r m") , s ") ( r c „// .,// 



ds ' / ds " E M' { ^ - MP , M")]) ( - 1 . r ( S ', M ') j 5 (r (*', /,') y ; r (s", //') , s") ( ( ^ ^ 

(66) 

In the Coulomb gauge, which we consistently use below, the most singular part of the 2x2 matrices C and D is the 
transverse part. Below we often only discuss this most singular part, i.e. the transverse contribution to S- la t, 



S intX x (*< t) = J2 A*V /"* ds' f ds" L 



int±±{ s \ A*'! S , //') 



(67) 



For later use we have incorporated the time interval (U , i) in which the interaction takes place into the notation for 
the action. The retarded Lagrangian L- lnt ±± = L- mt x±( s ' > M'i s " i A*") f° r a specific pair of paths only considering the 
transverse contribution is given by 



1 - ~ 'J 1 



£ £int_L_L = 6 (S - S ) - . 2 , 

ft « J (2?rft) 7 2tt 



rZq 



e -i«(y-«") exp 



-q-(r( S ', M ')-r( S ",M")) 



1 



q x r(s',fi') 



qxr( S " )M ") 



coth 



huj 

2k B T 



ImUn (q,w) 



q x r(s',fjf) 



dq 



do; 
2^ 



(27rft) 

qXp(V>") 



e -io;(s'-s") e i'q-[r(s',M')-r(s"./i")] 



(q,w) [l-2i<;(p (*"+,//'))] 



(68) 



In the second term ^-D_l_l can be replaced by either KeD±± or Im.D_i__i_ as long as we remember L lnt is retarded 
(s'>s"). 



G. Recapitulation of the major results 



Summing up, the linear response of the one particle density matrix to an external perturbation is given by Eq. 
which schematically can be written as 



(jw: 



cq 



This expression we approximately factorize in order to arrive at Eq. (63), or schematically 



6w= I (J)Z(V)Z- 



(69) 



(70) 



In the semi-classical limit the stimulus (-D)° q t is in the classical limit given by the driving term in the Boltzmann 
equation Eq. (^). In Eq. (^34j) the propagator (J)„ q t is represented as a double single-particle path integral with an 
action similar to the familiar Caldeira-Leggett influence functional. Schematically we write 



( J )a q t = / £>(path +) / £>(path -) exp 



+ Sint) 



(71) 



Here So is the usual action for non-interacting particles Eq. (^3l|). It does not couple the forward path + and the 
backward path — . Si nt contains the contribution from the gauge fields and the interactions and it couples the two 
paths. It is given by Eq. ([36|). The imaginary part of S- m t is governed by the fluctuations Eq. (|29|), while the real 
part is determined by the retarded gauge field propagator Eq. ( 28 ) . The distribution function also enters into the real 
part. 
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The Kubo formula 



stimulus 



-ffpert 



& 



propagation 



The path integral formalism 



stimulus 


& - 


D[/°] 





propagation 



J = J^V e is °l h 



The Boltzmann equation 

k 



stimulus 



D[/° 



propagation 



G- 1 = [r-V r + p- Vp]- 1 



response 



X = X H ] 



pert 



response 



X = JD[f°] 



response 



5f = G^Blf] 



FIG. 3. An overview of linear response theory for non-interacting particles in the Kubo formula, the path integral formalism 
and the Boltzmann equation. The linear response results from a stimulus combined with propagation. The figure illustrates 
where the distribution function f° appears in the three different cases. Note in particular the absence of f° in the propagation 
part of the path integral formalism and the Boltzmann equation. 



III. COMPARISON WITH THE KUBO FORMULA 



In equilibrium a standard way to perform linear response is to use the Kubo formula. Though the path integral 
formalism and the Kubo formula on a formal level are similar there are some important differences in the descriptions. 
This is most clearly illustrated by considering a non-interacting system as is done in Fig. [|, where we also compare 
to the classical Boltzmann equation. The linear response is in all three cases split into a stimulus that is propagated 
in order to give the response. In the Kubo formula the stimulus is the applied potential, while it in the path integral 
formalism is the induced excitation, which in the semi-classical limit equals the usual Boltzmann stimulus Eq. (|5C|). 
Despite that the stimulus is very different the Kubo formula and the path integral formulations have similar pictorial 
representations of the propagation: in both cases there is a forward and backward propagating electron line. However 
in the Kubo formula a liruipiathematically represents a Green's function and the propagator in the non-interacting case 
is the Lindhard functiontj, while in the non- interacting path integral formalism a line is a conventional free particle 
Feynman path integral. In the semi-classical limit the double Feynman path integral reduces to the classical path 
integral propagator for the Boltzmann equatiorO. The difference between the Kubo formula and the path integral 
formalism is high lighted by the observation that the (Fermi) distribution function in the Kubo formula is contained 
entirely within the propagator, while it in the path integral formalism is the stimulus and not the propagation that 
depends on the distribution function for a non-interacting system. Notice also that in the path integral formulation no 
statements are made about equilibrium. The fact that for a non-interacting system the propagation is independent of 
the actual distribution function is well known from the Landauer-Biittiker scattering approach to transport through 
a non-interacting part of the system. 

Within the Kubo formalism it is very natural to treat the interactions perturbatively. The effect of the interaction 
is traditionally divided into self-energy contributions, which alters the propagation properties of a single line or 
Green's function, and vertex corrections that connects the two lines. It is possible to make the same distinction in a 
perturbative treatment of the interactions in the path integral formulation, and we will also use the terms self-energy 
and vertex corrections to distinguish between interactions that only involves one path and those that connect the 
forward and backward propagating paths. For an interacting system the propagator in the path integral formalism 
thus also depends on the actual distribution function. This is well known in the classical limit of the Boltzmann 
equation, since the propagator then depends on both the collision integral and the Landau interaction function, but 
still the stimulus has already gathered the most trivial dependence upon the distribution function. 
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IV. THE SEMI-CLASSICAL LIMIT AND THE SEPARATION BETWEEN THE TWO PATHS 



Traditionally the semi-classical limit of a double path integral is the limit, where the important contributions are 
those where the forward and backward paths are almost identical. In a discussion of the semi-classical limit it is thus 
important to know how far apart on average the forward and backward paths are. This information is contained in the 
moments of the difference variable Ar (s). For a non-interacting system or more generally a system with instantaneous 
interactions the standard deviation of Ar can be estimated from the Heisenberg uncertainty relations between Ar 
and its conjugate variable, the average momentum P. That P indeed is the conjugate variable can be read off from 
the relations 

J^-So (Rf , P f , t; Ri, P i; 0) = -Ar (t) (72) 

^-S (R f , P f , t; Ri, Pi, 0) = Ar (0) , (73) 

where So is the action for a non- interacting system, Eq. (|6l|). We are mainly interested in Fermi liquids, so we prefer 
to use a polar coordinate system representation of the average momenta 

p w -'<•>(;;$)■ P4> 

where the angle along the Fermi surface <f> is measured relative to an arbitrarily fixed x-axis. The conjugate variable 
to the size of the momentum is 

■^-So (Rf , Pf , t; Ri, Pi, 0) = — IPf • Ar (t) . (75) 

From the Heisenberg uncertainty relations we then expect the standard deviation iP • Ar of the difference along the 
path to be given by h divided by the uncertainty in the size of the momentum, i.e. -pP • Ar rs Kv-g {kg,T) 1 . This 
argument is carried through in more detail in appendix ^ and the coefficient to hv-p {k^T) 1 is also shown to be 1/ y/3. 
The conjugate variable to the angle <j> is 

■^-S (Rf , Pf , t; Ri, Pi, 0) - -P £ x Ar (i) . (76) 

Remembering that the relevant momentum Pf approximately is the Fermi-momentum we will denote the quantity 
Pp 1 Pf x Ar by the width. The Heisenberg uncertainty relation then says that the standard deviation of the width is 
approximately the Fermi wavelength Ap divided by the angular uncertainty in the propagation direction. For a well 
collimated beam the width is thus big, while it is small (of the order of Ap) when measuring the current response to 
a homogeneous field, since in this case the angular uncertainty is as big as possible (of the order 1). 

The separation between the paths is not in itself directly measurable, but it tends to show up indirectly in experi- 
ments probing the semi-classical regime. It marks e.g. the cross-over between the long wave-length regime, where the 
vertex corrections are comparable to and partly cancels the self-energy and the short wave-length regime where the 
vertex corrections play a minor role. For instance in the calculation of the inverse transport time due to impurity 
scattering in the famous factor (1 — cos 6) the scattering angle 8 can be interpreted as j-Q^f, where q is the scattering 
vector and, as we saw above, Af is the width. 

For the composite fermions, where the force is transmitted via a magnetic field and thus a flux through an area, 
we expect the width to mark the regime where the vertex corrections sets in. We can thus use the width as an 
effective cut-off in the conventional calculations of the effective mass. The above considerations regarding the width 
are only true for a system with a time local action. For composite fermions this is not the case, since the important 
transverse fluctuations have a characteristic frequency around giq 2 , see Eq. (|33"1). We thus have to take the lowest 
order correction due to the terms in the action that are not time-local into account. We do this perturbatively. 

To this end let us define the average width squared (h/q^) 2 at time s as 



h^ 2 



PF 2 [P( s )xAr( s )f| [M , (77) 



where the average ^X^r 0)t i is defined as 
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W[<M] 



J ^ J dRf J ^ J rfRj ft (P f , R f ) D (Pj, RQ (X exp [ * S (R f , P f , t; R i; P i; 0)] )) [0 ,, 
/ A J dR f / A J dR, n (P f , R £ ) D (Pi, Rf) (exp [j-S (R f , P f , t; R h Pi, 0) 



(78) 



no,t] 



Here (X))[o.t] is a short hand notation for the path integral over paths starting at time and ending at time t with the 
Wigner function variables as boundary conditions. (-^)[o,tl is defined such as to include the action Si n t(i, 0) arising 
from the interactions within the interval [0, t] 



[o,t] 



[ VR [ VP [ VAp [ VAre* Sint{tfl) X 

J[0,t] J[0,t] J[0,t] J[0,t] 



(79) 



When no confusion is possible we will the index [0,t]. In the semi-classical limit the stimulus D (Pi,R;) is given by 
Eq. (p0|). For a current measurement the response 7£(Pf,Rf) = — evffi -1 = — e^-pffi -1 . Q is the volume of the 
sample. The time t is the time-interval during which the paths evolve. Below we will show that in the long time limit 
the cut-off is independent of t and s. We only consider this limit below and for our purpose t can be taken to infinity 
or rather be put equal to a scattering time. Notice that according to Eq. (164) the last factor in the denominator is 



exp 



— So (t, Rf , Pf , t; Ri, Pi, 0) 

n 



= (J(Rf,Pf,t;Ri,Pi,0)) at . 



(80) 



[0,t] 



It could easily be argued that it is better to use the absolute value of both the response 1Z (Pf , Rf ) and the stimulus 
D (Pi,Ri) in the definition of §X§, Eq. (|7q) , since the response might be the sum of a positive contribution from 
paths originating in one part of the stimulus phase space and a negative contribution from other parts of the phase 
space. The width squared of the positively contributing paths taken as a group is then canceled by the width of the 
negatively contributing paths taken as a group, the overall result might even be a negative width squared. This does 
not seem to be terribly important in the examples discussed below so for simplicity we stick with the above definition. 
We next use the conjugate variable relationship Eq. (f76|) to perform the rewriting 



pp 2 [P (s) X Ar (s)] 2 e tSa[R fl P f) t;R l ,P, 1 0] 



d 



So[R(s),P(s),s;Ri,Pi,0]. 



d 



s £So[Rf,P t ,t;R(s),P(»,s] 



ID 



dej) (s) dej) (s) 

inside the definition of the width squared Eq. (|77|). The action due to the interactions S[ n t is split into three parts 

rt 



(t,Q) 



V M y f ds' [ s d S "L int 

, „« = + Jo Jo 

v v'fi" [ d S ' r d S "L inu 

, „„ l Js Jo 



Sint (S, 0) + S int (t, S 



(82) 
(83) 



where the Lagrangian L ult only considering the transverse contribution is given by Eq. ( |68| ) . In total 

2 



-^1 =^ 



dP { 

{2irky 



dR 



J Eu',u»=± ds ' So ds " L 



dpi 



dRi^(Pf,R f )D(Pi,Ri) 



e |So[R(s),P(s),s;R 1 ,P,,0] 



dP{s) 
(2tt?i) 2 



dR(s) 



\\ dcj)(s) 
where the normalization is contained inside the new constant M 



dcj>( s ) 



f-S [R f ,Pf,t;R(s),P( s ),s] 



l[0,a]H[s,t\ 



dPf 

(2tt?i) 2 



c?R f 



dPi 



dRi K (Pf , Rf ) D (Pi, Ri) (exp 



-So(R f) Pf,t;Ri,Pi,0) 



))[o,t] 



(84) 



(85) 



The contribution due to the interaction is treated perturbatively; in particular we perform the expansion in 
the part that contains the scattering events that begin before the observation time s and finish later than s: 
E^= ± MV' />' I S ds"L int . 
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A. Zeroth order or the time-local limit revisited 



To lowest order or in the so-called time-local approximation, the last term in Eq. (|83| ) vanishes and the width 
squared is 

l)^ A ^/i^/ dR '/w/' iR ' R(p '' R ' )D(p " R ' ) /iw/ dR(s) 

(J [R («) , P (s) , s; R, P i; 0]) M (J [Rf , P f , t; R (s) , P (s) , s})^ . (86) 



X 



Let us study two distinct examples. The first is the case of homogeneous driving, where D (P;,R;) = D (Pi). We 
furthermore assume the system is rotation-invariant such that 

J dRi (J [R (s) , P (s) , a; R, P i; 0]) Oit (87) 

only depends on the size of the momenta Pi, P (s) and the relative momentum angle 4>(s) — fa. In this case it is 
nice to expand the angular dependence of the stimulus on the eigenfunctions exp (ivfa) for the angular momentum 
operator in the z-direction 

D(P0 =^exp(W> i )D I ,(P 1 ). (88) 

The width squared from the f'th mode is then i^ 2 Ap. In a path integral calculation of the impurity scattering integral 
in the Boltzmann equation the well known factor 1 — cos(yd) factor originates from an average of 2 sin 2 (g^ g^j j > so 

as expected it is the width z^Af which sets the boundary between the small angle scattering regime, where the vertex 
corrections given by cos(i^0)) are important and a large scattering regime, where the vertex correction average out. 
Notice that for an elastic scattering system this averaging has to be performed via the boundary conditions, while in 
general an interacting system is self-averaging. 

As another example we consider the gedanken experiment typically considered in connection with the interpretation 
of a single-particle Greens function. A particle is fed into the system at the Fermi surface and in a direction parallel to 
the x-axis. In order to make it a bit more realistic we give the direction a finite spread A; i.e. we take 1Z cx exp (^■^zfa 2 ^ ■ 
A time t later we measure the probability the particle is still moving along the x-axis with a directional spread of Af 
i.e. we take D cx exp ( fa? J . The width squared is in this case Ap (^-^ + -^z^J ■ In particular the width diverges as 
the angular spread to the power —1. 



B. The first order non time-local correction 



Let us move on to discuss the first order correction to the width squared due to the non time-local part of the 



interaction Y l ^', ll "=± /"V f a ds' f* ds"j- L h 



. i I dP { 
KJ ~(2nhy 



dPi 



A [-^-) =A F A^ / -jP^-z I JR { I / dRift(Pf,R f )D(Pi,R) / / dR(s) 



{2-Khy 



dP (s) 



£ MV" J* ds' jT ds"^L int (s', //; s", //') 



9 e iSo[R(s),P(s),s;Rd,Pi,0]_^_ fi iSo[Rf,P£,t;R(s),PCs),s] 



d<l>{s) 



d<f>(s) 



[0,s]//[5,t] 

(89) 



We have here neglected the contribution from the expansion of the normalization J\f, since it turns out to cancel a 
term in the expansion of Eq. (|8^) we anyway neglect in our discussion below. Let us here specialize to the contribution 
from the long wave-length fluctuations, where q • Ar -C 1, i.e. 



~L int ( s ', M W) 



If dq f duj , / hid \ . 

— / o / — coth haD ± j_ (q, 

hj {2irhf J 2tt \2k B Tj X± V4 ' 

Xe -ia)( s '- S ") e i-q-[R( s ')-R( s ")]+iiq-[ /i 'Ar( s ')- M "Ar( s ")] 



-qxR(s') 

q 



-qxR(s") 

q 



(90) 



15 



To lowest order in the interaction the extra contribution from L lnt to the width squared is in the long wave-length 
limit 



A 



/^/ dR '/^/' iR ' R(p '' R ' ,D(p '- R - ) s/S?/ dRW 
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To lowest order we will neglect the effect of Si n t (s, 0) and S; n t (t, s) inside the path integrals - 



and 



fat] 



(91) 

respectively. 



The product Ar (s') Ar (s") we will replace by its average value. The most important contribution comes from the 
longitudinal part of Ar (s') Ar (s"). Furthermore it turns out that the relevant time-scale in the above integral 
% (vpq)~ for the important fluctuations is smaller than the correlation time of Ar (s') Ar (s"), which is of the order 
Ti (fceT) -1 h ^p* ~ fi (/cbT) -1 , so we replace Ar (s 1 ) Ar (s") by its equal time longitudinal average value | ^ ^ , see 



appendix |^. Within these approximations the correction to the width squared is 
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For the composite fermions the spectral function -ImDn (q, ui) 
is approximately 
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Eq. (33), so the width squared 
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Here 77 = 1 corresponds to our prime example, the Coulomb interaction 
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(94) 



We expect the Fermi energy ^fPf to be of the same order of magnitude as the energy scale set by the interaction 

2 .„ / \ 3/2 3/4 

i.e. the average Coulomb energy £c = £^V™- The factor ( "J^f ) § (~) is thus of the order 1 and the width is 

approximately given by the geometric mean of the Fermi wave-length and the de Broglie wave-length. 

In perturbative calculations of the self-energy the cut-off is often formulated in terms of a frequency cut-off. The 
frequency cut-off corresponding to Eq. (B4J) is 



9 /tt\5/4 / 2e c 



V v f pf 



(95) 



We want to emphasize that this transverse cut-off is considerably larger than the longitudinal cut-off, which in 
appendix |^ is found to be \ // 3h 1 k B T. In the next section we show that in case of Coulomb interaction the effective 
mass depends logarithmically on w cut , so the temperature dependence of w cu t and thus the inverse width squared is 
VT I logT| 2 , which tends much slower to zero than the more conventional linear temperature scalingEl. 



V. THE EFFECTIVE MASS 

In the last section we calculated the width between the forward and backward path Eq. (^), and we hinted 
at how this width might play the role of an effective cut-off in the calculation of the effective mass. In this section we 
elaborate a bit more on this connection. 
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Experimentally the effective mass of a particle is determined by simultaneously measuring the momentum and the 
velocity. In a Fermi liquid the momentum is typically the Fermi momentum. The velocity measurement typically 
involves measuring how far the particle has propagated in a known time-span t. This is a little bit at odds with 
the semi-classical limit, where the external fields are considered to be smooth, so it is hard to define, yet measure, 
distances shorter than a transport mean free path length. This is necessary in a mass measurement, since the time 
t must be comparable to or smaller than the transport time to make the path reasonably well defined. Having said 
this one might hope for an intermediate regime, where it is possible to at least put some limits on the effective mass, 
while still maintaining the semi-classical line of thought. . 

One attempt at this is done by Willett, West, and Pfeiffeilij, who uses a surface acoustic wave (SAW) to generate 
an external electric field with a wave-length comparable to the mean free path. They then apply an external magnetic 
field and looks for resonances between the cyclotron radius and the wave-length of the SAW field. In order to obsetve 
these resonances the cyclotron frequency cannot be orders of magnitude larger than frequency of the SAW wave&3. 
This leads to an upper estimate of the elapsed time and thus the effective mass. 

We will denote this sort of experiment as a semi-classical mass determination, since the forward and backward paths 
together perform the cyclotron motion. This is in contrast to experiments probing the single-particle energy spectrum, 
like Shubnikov-de Haas measurements. In such experiments according to the Bohr-Sommerfeld quantization scheme 
it is the phase of a single path that causes the resonance. 

Below we show that in the semi-classical experiments the width between the two paths or in other words the 
enclosed area acts as an effective cut-off on the contribution from the magnetic field fluctuations to the effective mass. 
In the Bohr-Sommerfeld resonance scheme the cut-off is supplied by the area enclosed by the path itselfO. This area 
is much bigger than the semi-classical area and the cut-off is correspondingly smaller. 

In order to extract the semi-classical mass from the path integral description it is in principle necessary to analyze 
the concrete experiment in detail. Luckily it turns out that in the semi-classical limit it is possible to come up with 
an approximately local relation between the velocity and the momentum. This implies the effective mass is almost 
independent of the concrete mass measuring experiment in the semi-classical limit. 

For a non-interacting system the derivative of the action with respect to the momentum yields the Hamilton 
equation for the velocity — 
momentum Ap (s) 



= 0. For the two-path action it is the derivative with respect to the difference 
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(96) 



that yields the relation between the velocity and the momentum, at least in a stationary phase approximation where 
dAp(s) = 0- O n average the stationary phase condition is exact, i.e. 



dS 



dAp (s) 



= 0. 



(97) 



Notice that after performing the integral over R (s) and Ar (s) inside the average as defined in Eq. ( |78| ) we have 
Ap (s + ) = Ap (s) — q and P (s + ) = P (s) — |/Ltq. The terms with \i — fx' in Eq. (|9^) i.e. those terms where the 
interaction connect either the forward or backward path with itself gives rise to that part of the renormalization of the 
effective mass that is attributed to the self-energy in the usual perturbative treatments. This renormalization mainly 
depends on the local properties of a single path and it is thus independent of which pair of paths is relevant in the given 
experiment. On the other hand the terms with \i ^ // in Eq. (pw, which are associated with the vertex corrections, 
depend on both the forward and the backward paths. The renormalization due to these terms thus depends on the 
exact experiment. For instance in an interference experiment the two paths are almost independent and the vertex 
contribution is expected to vanish due to destructive interference in the phases associated with each path individually. 
In the semi-classical limit, where the forward and backward paths typically are close to each other it is fortunately 
possible to derive a universal expression for the renormalization of the mass. Below it will be shown the interaction 
modes with a wave-length longer than the separation between the paths does not renormalize the mass, while the 
renormalization due to the interaction modes with a shorter wave-length is given by the usual self-energy contribution. 
A possible interpretation is that the particle cannot distinguish between the field from a long wave-length interaction 
and an external potential. 
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For simplicity we neglect the extra terms arising from the velocity dependence of the force, i.e. we treat the 
interaction as scalar. With this small neglect the average of the velocity, denoted vp = ORb, is 
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(98) 



where in the first term it is used that the average momentum <|P|) is the Fermi momentum pp = ppWp 1 vp (see 
appendix ^|). At least for a non- interacting system Ap (s) is of the order of the typical wave- vector of the external 
force field. Our interest is the semi-classical limit so we will assume Ap (s) • pf < k-gT and neglect its contribution, 
i.e. we put Ap (s) = 0. The average velocity vp is then 



VF = Ip F _ f ds > f _^S_ / ^ VFm ^p F D ±± (ci^)e^i R ( s ')- R ^e-^( s '- s ) 
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(99) 



We are mainly interested in the long wave-length limit, so we will neglect terms of order q 2 compared to terms of 
order q. We also assume that both |P(s)| and the corresponding velocity R(s) are distributed evenly around the 
Fermi surface. In this approximation the second term inside the parenthesis can be neglected. We assume the relative 
variation in the velocity around the Fermi velocity vp is small, i.e. the temperature is much smaller than Fermi 
energy. It is thus possible to approximate the path locally by a straight line i.e. R (s') — R (s) i=s vf (s' — s). In the 
long wave-length limit the transverse part of q • Ar is bigger than the longitudinal part. The averaging procedure is 
approximated by independently carrying out the longitudinal and the transverse average. In the longitudinal average 
we use the time-local approximation outlined in appendix |^, so the longitudinal average of ^ ± — -^w [P (s) — /iq] 
is approximately 
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where we have used vp = ^p-- In the transverse average we approximate 
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with the result that the effective mass m*, defined as pf = m*vp, is given by 
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(101) 



(102) 



Apart from the extra factor of [l — exp (— 5<? ( Tl9 2 )] this relation is identical to the usual relation derived from the 
self-energy. Thus the contribution from the short wave-length modes of the interaction to the renormalization of the 
mass is identical to the standard result, while the long wave-length modes do not contribute. The cross-over is at the 
width hq~^ and this is why we above have alluded to the width as an effective cut-off. 

When the long- wavelength expression for the propagator Eq. ( p3| ) in the case of Coulomb interaction, is 

inserted into Eq. ( 102 ) the integral is ultraviolet divergent. We use pf as an approximative upper cut-off and get the 
effective mass 



4 pi 1 , pi 
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(103) 
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For composite fermions the renormalization of the mass is quite strong so the effective mass 

2 



is approximately independent of the bare mass to. 
The prefactor 



r 3/2 



Pf 



Pf , Pf 
— log; 



7r 3 / 2 2m ec 
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(104) 



(105) 



where ec = £^V™ 1S the average Coulomb energy and we as a typical valuesEj have used n = 1.6 X 10 15 m~ 2 , e = 13eo, 
and to w 0.067 too, where eo and too are the vacuum values of the permittivity and the electron mass respectively. 
The argument of the logarithm is given by Eq. (B3|) so self-consistently the effective mass is 
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6.3, 



(106) 



where the estimate is valid at a temperature of ~ 0.13 K. The cut-off at this temperature happens at distances 
hq~^ 6Ap 40 nm. 

As expected the semi-classical effective mass Eq. (106) is smaller than the mass deduced from the activation energy 
and Shubnikov-de Haas measuremental3, actually about half as big. In the experiment by Willett et al.E3 the estimate 



Eq. (106) leads to a cyclotron frequency thatrSA the secondary resonance position is about twice as big as the used 
SAW frequency 10.7 GHz. A naive argumenita leads to the expectation that the secondary resonance can only be 
observed if the cyclotron frequency is much bigger than the SAW frequency. It is thus a matter of taste whether an 
'apparent inconsistency' between the geometric resonance experiments and the size of the effective mass arises. We 
here want to stress the word 'apparent' since Mirlin and W61fleE£l has shown that the secondary resonance can be 
observed even when the cyclotron frequency equals the SAW frequency. 



VI. CONCLUSION 



In this paper we have developed a single-particle path integral description of transport for composite fermions, 
thereby extending the Caldeira-Leggett formalism for a single-particle interacting with an environment. Here the 
environment is the Chern-Simons gauge field. Our approach is modeled over the work by Golubev and Zaikin on the 
electron gas without a magnetic field. As in the Caldeira-Leggett formalism there is not one, but two paths in the 
path integral, one propagating forward in time and one backward. A new feature compared to the Caldeira-Leggett 
formalism is that the distribution function enters into the real part of the single-particle action. The imaginary part 
of the action is determined by the fluctuation spectrum, just as in the Caldeira-Leggett formalism. 

In a perturbative description of the CF's the long wave-length magnetic field fluctuations give rise to an infrared 
divergence of the self-energy and the vertex corrections, separately. In a transport measurement these divergences 
are known to cancel. In our formalism the infrared divergence does not appear because the vertex and self-energy 
contributions are never split. Not only is the infrared divergence absent, but the renormalization in the infrared limit 
is negligible. This result is obtained in a natural way by working in a real space formulation, since the magnetic 
interaction acts via fluxes, i.e. areas. In the path integral formalism the size of the particle is defined as the separation 
between the forward and backward propagating paths. Only the magnetic field fluctuations with a wave-length shorter 
than this particle size contributes noticeably to the renormalization of the mass. In GaAs at T ~ 0.13 K we find the 
mass renormalization due to the magnetic field fluctuations to be roughly 6 times the bare mass. Furthermore we 
have shown that in the semi-classical limit the separation between the two paths, i.e. the size of the particle, scale as 
T _1 / 4 | logT| _1 . This implies the effective cut-off in frequency scale as roughly T 1 t 2 \ logT| 2 . As usual the mass scales 
approximately with the logarithm of the frequency cut-off. 

The linear response regime around equilibrium is the major focus of our work, but the formalism is still valid out of 
equilibrium. In the non-interacting limit, as shown in Fig. 0, the single-particle path integral technique possesses the 
feature that the effect of the non-equilibrium distribution function is gathered completely within the stimulus, while 
the propagation is identical to the propagation in equilibrium. This feature is shared with other out of equilibrium 
techniques i.e. the Landauer-Biittiker description of transport through a non-interacting part of the system and the 
Boltzmann equation. In fact the classical limit of the single-particle path integral formalism leads very naturally to 
the Boltzmann equation. The single-particle path integral technique is thus a promising technique for the description 
of interacting systems both in and out of equilibrium, especially it may prove useful for studying systems out of 
equilibrium not easily treated by other methods. 
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APPENDIX A: EQUATIONS OF MOTION 

In this appendix we are going to derive differential equations for the composite fermion contribution to the effective 
action S c s[a; t] and the density matrix p ac (A e ff — a) in the presence of the Chern-Simons field A c g— a by differentiating 
the composite fermion part of Eqs. ( |l9| ) and ([l6|). 

In order to differentiate quantities in a fixed gauge field background at time s it is necessary to let the Chern-Simons 
field and the Grassmann fields ijj and ijj be defined on Ct, where t > s. The effective action for both the Chern-Simons 
field and the Grassmann variables defined in Eq. (oil) is then 



where 



Sct#> a ; s l = Scs[a;*] - e J da' J dra (r, a')n + ScomFerl^, ip, a; s; t], (Af) 

ScowFa['ip^,a;s;t] = J da' J drip(r, a')ih-^tp(r, a') - J da' H [A cf[ - a]. (A2) 

The second term in Eq. ( |AlD can be gauged away for t > s' > s, since there is no Hamiltonian. The terms that 
depend on the composite fermion Grassmann variables are gathered in S'comFer['0j a ! s i A an d it is its contribution 
we are going to consider in this appendix. 

Notice the effective action Sqf in the definition of p ac , Eq. (|16|) , can be replaced by ScomFor, since the difference 
cancels between the numerator and denominator in Eq. Jig). 

1. The effective action 

The composite fermion contribution to the effective action S^a; s ] is 

S cfi [a; s] ~ S cs [a; t] + e J ds' J drAa Q (r, s')n = ~ logff », (A3) 

where Aao(s') = ao{s',+) — ao(s',bi) is the difference field and the short-hand notation 

((A)) = Jl>$J O^exp (~S ComFe S,Tp,a;s;t} \ X (A4) 

for the Grassmann integral has been introduced. The derivative of Eq. ( |A3[ ) is 

a s {Eq. (A3)} = -{l}- 1 (H [A eS (s)-a(s,+)}-H [A eS (s)-a(s,-)]}. (A5) 

Per definition of the path integral is the right-hand side equal to the expectation value of — |flo[j4e£f(s) — o(a, +)] 

— Hq[A c h(s) — a(s, — )] | in the presence of the contour gauge field A e g — a. It can be expressed in terms of the density 
matrix as 

d s Eq. (A3) = - J dr (eAa (r, s)p ac (A e fi - a; r, r; s) + Aa(r, s) • j a (^l cff - a; r, s)) , (A6) 

where the current in the presence of the gauge field A c r — a: j a (^4cff — a;r, s) is defined in Eq. (|2C|), Eq. (pl| ) then 
follows upon integration using the conventional assumption of adiabatic turning on of the gauge field. 
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2. The density matrix 



The density matrix in the presence of the Chern-Simons gauge field p ac is defined in Eq. (|l6|). At a time a < t the 
action Scf[iPi "0: a 'i s ] can be replaced by 5comFer[V ; j °j s ! £]■ ln terms of the above short-hand notation 



p ac (r,r' )S ) = fl)- 1 #(r' )S )^(r )S ))). 
When a is varied both the numerator and the denominator varies: 

iW sPac (r,r',s) = ll}-hhd s lHr', s )iP(r,s)} - {l}' 1 ^', s)^(r, s)}ihd s { 1] 



(A7) 



(A8) 

In the above subsection [A l| the derivative of the denominator ift<9 s {l)) was calculated. The contribution from the 
first term (the numerator) is 



We here used 



! {((/f L4 off (s) - a(s, +M(r', 8 )i>(r, s) - ^(r', aty(r, a)# L4 cff (a) - a(s, -)])}. 
#(r', a + As)^(r, s + As))) = f^(r', s)V>(r, a)), 



(A9) 
(A10) 

(All) 



since without a Ham iltonian the system does not evolve. The action <ScomFer is quadratic in the Grassmann variables ip 
and ip, so Eq. jA10| ) can be evaluated with the help of Wick's theorem. The Hartree disconnected diagram like pairings 
exactly cancels the contribution from the variation of the denominator in Eq. flAS|), fl))" 1 ((^(r', s)^(r, s)))iftd s ((l)), 
so paying due attention to the contour ordering 



ihd tPac (r, r', t) = /" dr" / dr'" (<5(r - r") - p ac [r, r"; t])H(r", r'"; f, +)p ac [r'", r'; t] 
- J dr" J dv"'p ac [r, r"; t]H(r", r'"; t, -) (6(r"' - r') - p ac [r'" , r'; i] 

This is Eq. @ written out. 
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APPENDIX B: LINEARIZATION OF THE DENSITY MATRIX DIFFERENTIAL EQUATION 



The stimulus in the linearized equation of motion for 6p ac (A e ff — a; r, r', t), Eq. (46), D = Di + D2 is split into two 
parts. The explicit form for Di, which is independent of Aa, is 
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(Bl) 



The explicit form of D2, which contains all the dependence upon Aa, is 
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APPENDIX C: THE LONGITUDINAL SEPARATION 

In this appendix we are interested in the distribution of the longitudinal separation -pP ■ Ar and its conjugate 
variable the size of the momentum P. For simplicity we will only consider the case of a scalar interaction. 

The first step is to rewrite the free particle action such that the longitudinal and transverse directions appear 
explicitly. 

To this end we use polar coordinates for the momentum 



P = 



psin — 
1 v 



(CI) 



— is the angle with respect to an arbitrarily defined rc-axis. p v , loosely speaking the arc-length along the Fermi 
surface, has been chosen as our variable instead of the angle in order to avoid a complicated Jacobian. In the 
scalar approximation the interaction is independent of p n , which therefore can be integrated out. The free particle 
action on the forward part of the contour is 



I dS (" P 



P 

2m 



The stationary phase approximation for p v reads 



p„ = p arctan — . 

x 



(C2) 



(C3) 



There are two solutions to this equation corresponding to the momentum and the velocity being parallel or anti- 
parallel. Here we only keep the parallel solution. The resulting stationary phase action on the forward contour 
is 



L 



where v — x 2 + y 2 is the velocity. A similar result holds on the backward path, but with the opposite overall 
sign. We are interested in discussing the separation between the two paths, so we change representation from the 
position and momentum along each path individually to the average R (s) = | [r (s, +) + r (s, — )] and the difference 

Ar (s) = r (s, +) — r (s, — ) . In a translation invariant system one has the choice between using R or R as the variables 



in the path integral. Here we choose R and we parameterize it using the velocity V(s) = R(.s) 
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and the polar angle 



(C4) 
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The unit vector n (s) points along the direction of the average propagation, n (s) can also be thought of as a unit 
vector on the Fermi surface, since above in Eq. (C3) we only used the stationary phase condition where the momentum 
is aligned with the velocity. 

Neglecting third order terms in the difference variables the free particle action is 



„ Jo 



p(a,n)v(s,fi) 



2m 
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, (C5) 



where k is a Lagrange multiplier vector that has been introduced in order to deal with the boundary condition in 
the transition from R to R. i|k|Wt r is a small number, when the feeding-in and the pick-up are not local in space. 
As an illustration assume that the initial and final position is distributed according to isotropic Gaussians centered 



around R? and Rf respectively. The sum of the variances is denoted by (-jx) 
is then replaced by 



The term k • (Rf — R;) in the action 
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after the integration with respect to Ri and Rf has been performed. The typical size of k is now given by the smallest 
of the numbers Rf ^ R .| and AA. The last number is in the semi-classical limit small. In the case of a periodic driving 
k is the wave-vector of the external field. Anyway in the semi-classical limit we will neglect the contribution from k. 



1. The interaction part: general formulation 



We first Taylor expand the interaction part of the action Si n t, the scalar part of Eq. (pq), to second order in the 
difference variable Ar. Furthermore, we use a time local approximation. That is we approximate the quadratic terms 
Ar(s)Ar(s') as Ar(s) 2 . Using this 
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F (R, a) a0 = {5 a0 - n an/3 ) F T [V («)] + n a n F h [V (t 



f* , „ dD 00 „ ^vf 1 (v>("+\,, Ap{s"+) 
J ~d\R\ ( V ( s >( s ~ 8 ) ' s ~ s ) 2-j [2 ~ W \ ' S ' ( ' 2 



= / ds 



00 



em 



(V(s)\s'-s\,s'-s) 



ds'- 



dC Q 



V(s)\s'-s\d\R 



(V{s)\s' - s\,s' - s) 



(C7) 



(C8) 
(C9) 



(CIO) 

(CII) 
(C12) 



As a further simplification we have above also approximated the average path by a straight line inside the kernels. In 
fact we will below treat both Ft and Fl as constants. 

In this approximation the action separate into a longitudinal and a transverse part. We define 



Ar (s) = Ar|| (s) n (s) + Arj_ (s) z x n (s) . 



(C13) 



In the Wigner variables basis of the propagator we use a polar coordinate representation of the initial and final 
momenta 
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Pi = ,pi 
Pf = P{ 



COS (pi 

sin 4>i 

COS (j)f 

sin (f>{ 



(C14) 
(C15) 



In the full propagator we now split out all the factors related to the transverse direction into the so-called transverse 
propagator 



Jt 



Pi;, / VAr ± ^5 



V(t) 



2-kTi 

IT' 



V(0) 



kf*da iF T Ar 1 ( s ) 2 +l Ar 1 ( s )P( s )£^ 



(C16) 



In terms of Jt the full propagator is after integrating out Ar(0) and Ar(t) 

J = j VP J VAp j VV J VVAr\\ 2ith5 [Pf - P (t)} 2nh5 [P - P (0)] J T 

xe i J o t d s [Ap( s )^( s )-ip( s )Ap( s )]-l f*d s [±Ar n ( s ) 2 F L +iA ril (s)(£;P(s)+hlV(s),P,s])} 



(C17) 



To lowest order the transverse propagator is independent of the dynamics of the longitudinal variables, so we will 
treat Jt as a constant and absorb it in the overall normalization of the path integral. The study of the remaining 
longitudinal propagator is the topic of the next subsection. 



2. The interaction part: the longitudinal contribution 



In the propagator Eq. ( |C17') we will neglect the dependence of h upon Ap. V (s) thus equals m 1 P(s). The 
propagation described by Eq. flC17) is then equivalent to the quantum propagation problem 



We neglect the momentum dependence of -FL- In equilibrium is 



h(P)ma- tanh 



2k B T 



(P-Pf) 



where the constant a in the high temperature limit is related to the fluctuations through 

fiF h k B T 



(C18) 



(C19) 



(C20) 



The last term in Eq. ( |C1SD is gauged away by the transformation <fr(P) ^— ^jfj and the propagation probler. 
Eq. (|C1SD is equivalent to 



M n 2 Fld\ 

2 — 7^0.6 = 



a 2 dP 2 



hF L d 
a dP 



tanh 



1 v-p 



2k B T 



(P-Pf) 



tanh 



1 Vp 



2k B T 



(P-Pf) 



(C21) 



The ground state of the quantum problem Eq. (C21) is the derivative of the Fermi distribution function and the 
ground state energy is zero. The typical decay time of the excited states is given by ^t. This can be identified as the 
correlation time scale of Ari or P with themselves. It is also the time scale it takes before the boundary conditions 

cannot be felt anymore. For the Coulomb interaction this time-scale is approximately h (k B T) 1 TigiPp (vppp) 1 . The 



ground-state wave-function of the original quantum propagation problem Eq. (CIS) is proportional to w — This 

implies the distribution of the longitudinal momentum or energy is given by the derivative of the Fermi distribution 
function, while the average of the width squared is 

-1 poo /-t d2/ \ 2 / a , \-l 



dp 



dfp 
dp 



dp 



h d 2 fp 
2 dp 2 



dfp 
dp 



i h 2 v 2 

3 (k B Tf 



(C22) 
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